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MOMENTS OF 9—NORMAL AND CONDITIONAL g—NORMAL 

DISTRIBUTIONS 

PAWEL J. SZABLOWSKI 


Abstract. We calculate moments and moment generating functions of two 
distributions: 

the so called q —Normal and the so called conditional q —Normal distribu¬ 
tions. These distributions generalize both Normal (q = 1), Wigner (q = 0, 
q— Normal) and Kesten-McKay (q = 0, conditional q— Normal) distributions. 
As a by product we get asymptotic properties of some expansions in modified 
Bessel functions. 


1. Introduction 

The purpose of this short note is to present exact forms of moments and moment 
generating functions (i.e. Laplace transforms) of four distributions with densities 
that are denoted by: fh{x\q), fN{x\q) and fq(x\a 1 b 1 q), f C N{x\y, p,q). In fact dis¬ 
tributions fh and /at are related to one another by the linear transformation of 
random variables having these distributions. Similarly distributions /q and /cjv 
are interrelated. The details will presented below. Two of the presented below dis¬ 
tributions are called respectively q —Normal (/at) and conditional q —Normal (/civ)- 
These distributions are the elements of the chain of attempts to generalize Normal 
distribution that exist in the literature. The first of the considered in this paper 
distributions (/at) was described first in the noncommutative probability context 
in [6], later in the classical probability context in [7]. Of course it was not the only 
one attempt to generalize Gaussian distribution. For others, different see e.g. [24] , 

m,m- 

Let us mention also that for particular values of the parameter q we get Gauss¬ 
ian (Normal) (q = 1), Wigner and generalized Kesten-McKay (q = 0) distribu¬ 
tions. Let us remind that distribution /at appears in many models of the so called 
q —oscillators that are considered in quantum physics (to mention only [4], [9], [3], 
El)- On the other hand Wigner or semicircle and Kesten-McKay distributions ap¬ 
pear as limiting distributions of certain combinatorial considerations and also in 
the context of random graphs, random matrices and large deviations. By the gen¬ 
eralized Kesten-McKay distribution we mean distribution that has density of the 
form CVa? — x 2 /Q 2 (x), where Q 2 (x) denotes quadratic polynomial that is positive 
on [— a , a] and C is some normalizing constant. Formal definition dating back to 
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papers of Kesten pTD] or McKay JT3] concerned very special form of quadratic poly¬ 
nomial Q. For more recent uses of Kesten-McKay distribution see e.g. m, mi, 

PS- 

The distributions that we are going to recall in this paper have appeared also in 
the context of stochastic processes allowing generalization of Wiener and Orstein- 
Uhlenbeck (see e.g. |JE]) processes and also in the context of quadratic harnesses 
(for review of the rich literature on this subject see [17l). 

The paper is organized as follows. In the next section we recall basic notation 
and basic properties of the analyzed in the paper distributions. In Section [3] we 
present our results. In Section [4] we collected longer proofs. 


2. Notation and basic notions 

To present these distributions we will use notation commonly used in the context 
of the so called q— series theory. Nice introductions to this theory can be found [X] 
or [8j. 

So q will be a parameter such that q £ (—1,1 > • For |g| < 1 the formulae will 
be explicit, while the case q = 1 will sometimes be understood as a limiting case. 

We set [0], = 0, [n] q = 1 + q + ... + q n ~\ [n],! = TI”=i \J] q , with [0] ? ! = 1, 



[n-k\ q \[k\ q \ 

0 


n > k > 0, 
otherwise. 


We will use also the so called q— Pochhammer symbol for n > 1 : (a; q ) = 

nj=o i 1 - ac t J ) > (ai,a 2 ,...,a k -q) n = n)=i With (a;q) 0 = 1. 

Often (a; q) n as well as (ai, 02 ,..., a,fc; q) n will be abbreviated to ( a) n and 
(ai, 02 ,..., cik) n , if the base will be q and if such abbreviation will not cause mis¬ 
understanding. 

It is easy to notice that for q £ (—1,1) we have: (q) n = (1 — q) n [n] q \ and that 
iq)n n > k > 0, 

otherwise. 


GL = 


(q) n -k(i)k 

0 


To support intuition let us notice that: [n^ = n, [rz] x ! = n\, [^] 1 = (^), (a; l) n = 


(1 — a) n and [n] 0 = 


if n > 1, 
if n = 0, 


W = 1, [*]„ = !.(«;'o) n = 


1 if n = 0 
1 — a if n > 1 


Let us denote for simplicity the following real subsets: 


J{q) 


[-2/vT^,2/yT^] if \q\<l 
K if q= 1 


( 2 . 1 ) 
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The four distributions that we are going to consider in the paper are defined by 
the densities: 


( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 


fh{x\q) = 2 ^)°° y/i - x 2 JJ((i + q k ) 2 _ 4g fc x 2 )/ [ _i i i ] (x), 

j=l 

f N (x\q) = 2 ft ((1 + q k ) 2 - (1 - q)x\ k ) I J(q) (x), 

27T\/4 — (1 — q)x 2 k _ Q 

c t \ u \ (l,ab)oo rr ((! + 9 fc ) 2 - ^q k x 2 ) , , 

f Q {x\a,b,q) = -— II - — j—;—r-/[-i,i](a:), 

2tt\/ 1 - a; 2 ^ w k {x\a,b,q) 

°° ("1 _ „2 k \ 

fcN (x\y, p, q ) = /a (a:|g) n w k (x,y\p, q) Ij{q) (x) ’ 


where we denoted: I a (x) = 
polynomials: 


1 

0 


if 

if 


x € A 
x A 


, and Wk and W k are the following 


( 2 . 6 ) 

W k (x, y\p, q) = { 1 - p 2 q 2k ) 2 - (1 - <?)pg fc (l + p 2 q 2k )xy + (1 - q)p 2 {x 2 + y 2 )q 2k , 

(2.7) t/;fe(a:|a, b , g) = (1 + aV fc )(l + &V fc ) - 2:r(a + b)q k ( 1 + a6 9 2fc ) + Ax 2 abq 2k . 
k = 0 , 1 , 2 ,... . 

Notice that Vfc > 0 : w/c (cc|a, b, q) = ui 0 (a;|ag fe , bq k , l) , W k (x,y\p,q) = Wo(x, y\pq k , q) 
and that Wk (x,y\0,q) = 1. 

Parameters characterizing these distributions (other than q) have the following 
ranges: y £ J(q), \p\ < 1, |a| , |6| < 1. 

These densities are defined for |g| < 1 with possibility to extend this range to 
q £ (—1,1] for densities /jv and Jcn but the cases q = 1 will be understood as limit 
cases. Thus important special cases can be summed up as follows: 


fh(x | 0 ) 
fN(x |1) 
fcN{x\y,p,0) 


7 Vl - x 2 /[_i,i] (x), /w( x| 0) = ^-y/4-a: 2 / [ _ 2i2] (a;), 


- 7 == exp (-x 2 /2) , /q(x|oi, b, 0) 

v Z7T 


2(1 — ab) \/l — x 2 
nwo(x\a, b , 1) 


(1 — p 2 )\/4 — a: 2 
27rW 0 (a;, y|p, 1) ’ 


/cjv (z|y,p, 1) 


1 

VMl 



(a-’ ~ ffl) 2 

2(1 P 2 ) 


It is known (see e.g. 0 but also detailed review [21]) that these distributions make 
the following families of polynomials orthogonal. These families will be defined 
through their 3-term recurrences: 


(2.8) h n +i(x\q) = 2xh n (x\q) - (1 - q n )h n -i(x\q), 

(2.9) H n+ i(x\q) = xH n {x\q) - [n\ q \H n - 1 (x\q), 

(2.10) 

Q n+ i(x\a,b,q) = (2x - (a + b)q n )Q n (x\a,b,q) - (1 - q n )( 1 - a&q'" _1 )(5„_i(x|a, b, q), 
(2.11) P n+1 (x\y,p,q) = (x- pyq n )P n (x\y, p,q) - (1 - p 2 q n ~ 1 )[n] q P n _ 1 (x\y, p, q), 

with h-i(x\q) = H-\(x\q) = Q-i(x\a,b,q) = P-i(x\y, p,q) =0, h 0 (x\q) = H 0 (x\q) 

= Qo{x\a,b, q) = Po(x\y, p, q) = 1. Polynomials h n and H n are called q— Hermite, 
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(more precisely continuous q— Hermite), while Q n and P n Al-Salam-Chihara poly¬ 
nomials. 

It is also known that the orthogonal relations have the following form: 


(2-12) J H m (x\q)H n (x\q)f N (x\q)dx = j jjj , ^ n = m ’ 

(2.13) J ^h n (x\q)h m (x\q)f h (x\q)dx = j ^ ^ m = n ’ 

(2.14) J ^Q n (x\a,b,q)Q m {x\a,b,q)f Q (x\a,b,q)dx = j ^ ^ ™ 7 ” 


(2.15) 

/ -Pm(a:|y,P,g)-Pn(a:|y,P,g)/civ(a;|y,P,g)da; = 

Jj{q) 

There are interesting special cases: h n (x |0) = U n (x), H n (x |0) = U n (x/ 2), 
ff„(a;|l) = -ff„(a;),Q n (a;|a,&,0) = t/„(x) - (a + b)U n -i(x) + abU n - 1 (x), P n (x\y,p,0) 
= U n (x/2) - pyU n -i(x/2) + p 2 U n -i{x/2), P n (x\y,p,l) = (1 - p 2 ) n l 2 H n (fx - 
py)/V± — p 2 ), where we denoted by U n n— th Chebyshev polynomial of the sec¬ 
ond kind and H n (x) denotes ordinary Hermite polynomial (so called probabilistic) 
i.e. monic orthogonal with respect to measure with the density: exp(— x 2 / 2) /\f2it. 
All the above mentioned facts can be found in [8] but also in more detail in [5], 

[?]. na. Ca¬ 

in the sequel we will need the following two facts: 

(2.16) J h n (x\q)f Q (x\a,b,q)dx = S n (a,b\q), 


0 if n m 

(p 2 )„[ n L ! */ n = m 


(2.17) [ H n (x\q)f C N(x\y,p,q)dx = p n H n (y\q), 

JS(q) 

where S n (a,b\q) = X)”=o The first °f them is shown in [20], the second 

in [5]. Notice that S n (a, b\l) = (a + b) n , S n (a, 6|0) = (a n+1 — b n+1 )/(a — b) if a ^ b 
and S n (a,a\0) = (n+ 1 )a n . 


3. Moments 

We will need the following and the following expansion: 


Ln/2J 

(3.1) (l-g)"/V= Yl ((") 

k =0 ' 7 

that can be easily obtained from the relationship: 

xU n {x/2) = U n+ i(x/2) + U n - 1 (x/2). 

This expansion can easily be modified to obtain the following ones: 


n 

k- 1 


)U n - 2 k(x\/l - q/2), 


(3.2) 


L ^/ 2 J / \ / \ 

2 ”*” = E((0 -U-! “)' 

k =0 v 7 V 7 

L n /2j / \ / \ 

xn = £ ( Q ~ (k - i) ) c/ "- 2fc ( x / 2 )- 


(3.3) 
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Let us remark that (n — 2k + 1) ("£ 1 ) /(n + 1) = (£) — ( fc " 1 ). 

Basing on these expansions we are able to formulate the following proposition 
giving expansions of x n in the series of q— Hermite polynomials. 

Proposition 1. Let us denote c m , n {q) = (- 1 ) J g j(j+1)/2 ((Jl J .)-( m ^._ 1 ))[ n_2 ” l+J ] 

defined for n > 1, m < |_n/2j . Then 


Ln/2j 


\n/2\ 


^ ' (1 q) “ e rn ^ n H n —2m{x\qf x ^ ( c rn ^ n h n —2m (*£k)• 


m=0 


m =0 


Proof. We use AMD and (E2) and the two following expansions first of which was 
shown in jTB] (4.2) and the other is its obvious modification: 


L«/2J 

U n (x^l—q/2) = (- 1 )' t 1 - «)"/ 2 -V W+1)/2 

i=o 

Ln/2j 

t^(*) = (- 1 )V ' 0+1)/2 

i=o 


n-j 

L J 


H n _ 2 j (x\q) 


J 9 


n — j 

. j . 


(x|g). 


□ 


Proposition 2. Let A' ~ /jv ttien i) Vn > 1 : 

(3 - 4) _ 

0 if n is odd 

■Ei=o(-l) fe '? m (2fc + l)m i ) if if n = 2j 


(1 - q ) n/2 [ y n f H {y\q)dx = 


2j+l 


rr—Z — 00 _ 

(3.5) <p N (t\q) =Eexp(t,X) = -^^(-1)V ^ )(2 k + l)J 2fe+ i(2t/- v /l - g), 

fc=o 

where Ik (t) is t/ie modified Bessel function of the first kind, 
ii) Let X ~ fh then 

(3.6) 

0 if n is odd 

4JWq^ELo(-l) fc 9 (T) (2fc + l)(^ + fe 1 ) i/ n = 2i ’ 


a" A (a|?)da; = 


2 fe-(-1 

V/^k) = -Eexp(tX) = y^(-l)V = )(2fc + l)/ 2fc +i(t). 


Proof. Is shifted to Section 2) 
Remark 1. Setting q = 0 we jet 

/ J 2n /ff (y|0)da- = 


k—0 


□ 


2 n 

n — 1 


2?r 


n + 1 V n 


Le. n—Catalan number. Setting q = 1 we get Vn > 1 : 

1 


(2n—1)!!= lim 


?->i- (2n + 1)(1 - j)* 


^ ( -i)fe 5 m(2fc+i) 


k—0 


2n + 1 
n — k 


As far as the moments of /cat and /q are concerned we the following result. 
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Proposition 3. Let X ~ /cjv with parameters y,p,q then 


i) 


\n/2\ 


EX n = £ (1 - 9r m p"- 2m tf,- 2m (yk)c m ,n(g), 


m =0 


VcNitiViPiV) = Eexp(tX) = g ^ ^ — p k H k (y\q) 


t ^ [k)J 

fc=0 L J 9 


x E(-ir [t+jV |* + f 3 + V »v.wVh). 

j=0 


«) Let X ~ /q with parameters a, 6, g £/ien 


Ln/2j 


EX n = 


^ ^ (o, 6|^), 


2=0 


= j g; w<^+ + 0 grtw)/* W ,(t). 


^ fc=0 ^ k j= 0 


Proof, is shifted to Section [I] 

As a corollary we get the following relationship. 

Corollary 1. i) lim^- c m ,„(?)/( 1 - q) m = 

1 


□ 


lim , 

9^1- (1 - 9) r 


£(-!)V 


2=0 


2„2(2+l)/2r/ ^ 

m - j 

n\ 


n 

m — j — 1 


n — 2 m + j 

j 


J 9 


2 n m\(n — 2m)! ’ 




_ OO _ 

lim -—^£(-l) fe g( * ) ( 2 k + l)/ 2fc+ i(2t/y/l - g) = exp(-t 2 /2), 

<j—>1“ t z —' 


fc =0 


mj 


r (1 — ?) fe ^ 2 k tt / i \ v~v 1 u [^ + + 2j + 1) 

llm ^^E m , P fc gfe(2/|g)£(-l) J - 2 - 




2=0 


[iU 


X q 2 ( 2 + 1 )/ 2 / 2j+fc+1 ( 2 i/ x /i _ g ) = exp(tpy + (1 - p 2 )t 2 /2). 
Proof, i) We apply the following well known expansion 

L"/2J 


= E 


—' 2 n m\(n — 2m)! 


H n —2m (#) 


m —0 


and the fact that f_ QO H n (x)fcN(x\y , p, l)<ir = p n H n (y) obtaining n—th moment 
of the fcN(x\y, p, 1) distribution: 

L n /2J , 


/ °° _ z 

x n f C N{x\y,p, l)dx = £ 

■°° m=0 


2 n m!(n — 2m)! 


H n — 2m(y)- 
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Now applying uniqueness of the expansion in orthogonal polynomials and assertion 
ii) we deduce our limit. 


ii) We have J^ exp (tx) exp(— (x — py) 2 /(2{\ — p 2 )))dx/^/2n(l — p 2 = exp (tpy + 


(l-p 2 )i 2 /2). 


□ 


4. Proofs 


Proof of Proposition [H i) We use the following expansion following expansion and 
its obvious modification: 

r - oo 

(4.1 )f H (x\q) = ^4 ~ (1 - q)x 2 J2 (-l) m 9 r2+1) C/2 m (xVT^/2)), 

ra=0 

9 j -i 2 < ^ > 

(4.2) f h (x\q) = ~ X ]T(- ir q ( m n U 2 m (x), 

7 r ' 

m=0 


given in If5] [Lemma 2, iv]. Secondly we apply (13.11) and the fact that polynomials 
U n {x/ 2) are orthogonal with respect to Wigner distribution, 
ii) We have: f^ 2 ^l-exp(yt)f H (y\q)dy = 

E OO 

7=0 


(24+ !)&') 


■‘3=0 (2j + l)!(l —q) 

= Er=o(-l) fc 7 (fctl) (2 k + 1) TZk 773 


pi 


jk— 0 V 

= Ek=o( —2 ) (2fc + 1) Em=0 m!( 2 fc+m+l)!(l-q) fc +” 
that I a (t) = E™=o Jrff+a+i) • 

To show (GDI we have: exp(xt)f h (x\q)dx = E^lo ' ( 2 j+i)I 4 j" Efc=o(~ 1 ) fc g^ 2 ^ ( 2 ^+ 

!)(?+*)=Er=o(-i)^ (tll) (2Hi) e~* (J - fc)! j+ fc+1)! 42 = Ez 0 (-^ m ^+ 


■J=fc (l-fe)!0'+fe+l)!(l-9) 3 


Now it is enough to recall 


i)E“= 


0 ra!(2fc+ra+l)!4 fe + n 


□ 


Proof of Proposition^ i), ii) We use the following three facts: One is m, the 
second (12.171 l2T6l) . the third formulae (14.11) and (14.21) . Using them we have: 

(1 -q) n,2 EX n = ^ T Efc=o 2J («-2fc+l)("fe 1 ) I%'f^ i U n-2k{.Xy/T~ql2)fcN{x\y,p,q)dx 

= 7^Ei”= / o J («-2fc+l)("+ 1 )x 

Ejn/'J-fc ( _ 1}J (1 _ qT l2-k- 3qj{ j+l)l2^-2k-j^ pn -2k-2 ]Hn _ 2k _ 2] (y | g) 

= ^rElto J (i - g) n/2 - m p”- 2m i?„- 2m ( 2 /|g)x * 

E7=o(-l)V°- +1)/2 (n - 2 m + 2 j + 1)(^) 7 +J ] 

2/> P, <l) = E7=0 £ Em=0 t 1 - g)- m p”- 2m .ff„_ 2ro (y|g)c mi „(<?). 

We have further: 


PcnZ y, p, q) = E“= 0 i 2m /(i - 7) m E7=2 


jrn—0 L / \ ± V) Z^ n =2m 
4.2m //i „\m Y^°° £ fe „/c i 


Hn—2m (y\Q)Cm,r> 


~ Em^O ^ 2m /(1 ^) m S/c=0 (fc+2m)! P Hk{y\q)c rri ,k+2m((l) ~ 

Efe =0 Cy Hk{y\q) Em=o i ■ i /-j... i! 

=Er. n tygy*) e”-o n-,)-‘w„ + i)t E J m -.(- i ) , y 0+1)/2 ( t +2j+ i )( t t;? 1 ) py], 

= E7=o wiP k H k (y\q) E°l 0 (-1 V oo 


_ J KJ , ^ Y^OO _ 

J=0V •‘■7 (l-9) J [i]J 9 2^m=j (1 -g)" 


(ra—j)!(/c+m+j + l)! 
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_ Y^°° t k „k U (nAnW^ 00 < i M' ' „i(i+ 11/2 V^oo _ t^_ _ 

2-^k—0 [k\ q \r kyyiV) 2-*/j=Q\ ) [j] g !(l—gR " 2^n =0 (l—q) n n\(k+2j+l+n)\ 

= TZo [ife P k H k {y\q) E^lo(- 1 ) J ** [k (i-ivmj +1) Q jU+1)/2 
x(l - q y+ 9k / 2+1 /H-^- k -H 2j+k+1 (2t/VT^) 

= ^ Er=o fefeV^(yk) E” o(-i) j [k+i]q '^ 2j+1) ^-o - +1)/2 / 2j+fc+ 1 ( 2 t/vr^g) 

ii) The proof of the first formula is analogous. 

Further we have: 

a, 6, g) = £~ = o Ej=o 2J c jtn S n - 2 j(a, b\q) 

= E”=0 fe E“=2m 

= Em=0 4™" En=2m 2 fc (fc+2m)! ^ fc ( a > &l9) c m,fc+2m 
= Er=o Em=o ^c m>fc+2m /(fc + 2m)! 


E oo t fc S fc (a,%) v-^°° 

/c—0 2^ Z-/m=0 4 71 


__ Y^°° t k S k (a,b\q ) Y^oo [ Hi t 21 (q) k+1 {k+2j+l) j(j+ 1)/2 v^ 00 

— Z^fc=0 2 k (q) k Z_y,=0V V 2 2 JfaU ^ Z^m= 


■ Er= 0 (-i)v w+i)/2 (fc+2j+i)( fe+ ^E +i ) m 


j 

u.2m — 2j 


2 2 H<3h 


m=j 4 m -^ (ra— < 7 ')!(fe+m+j+l)! 

t 2n 


_ V^oo t fc S fc (a,b|g) y^oo / i J (g) fe+J (fc+ 2 j+l) 7 ( 7 + 1)/2 y^oo _ 

2 fc (g) fc Z-^j=0\ - 1 / 2‘ 2j (q)j V Z-^/n—O 4 n (n)l(k+n+2j+l)l 

_ Y^°° * fc Sfc(a,%) t i M' * 23 (g)fc+3(fc+ 2 j + l) „jf, + i 1/2 2 2J+fc+1 r 

— 2^fc=0 2 fe ( 9 ) fe Ej=of -*-1 2 2 J(g) 3 - 9 t 2j+fc+l -^j+fc+l It J 

— 2 W°° g fc (a,b|g) Y^o° / IY7 (g)fc+J(fc+2j+l) ? -(. ? - +1 )/2 r . , , , , fyV n 

— t 2^fc=o ( q) k 2_/j=ov i l ( g ) 3 - y r 2 j+fc+urj- i-i 
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